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We study the interior of distorted static axisymmetric blacfc holes. We obtain a general interior 
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I. INTRODUCTION 

The uniqueness theorem proved by Israel [l| tells us 
that the only static vacuum black hole solution of the 
Einstein equations in an asymptotically flat spacetime is 
the Schwarzschild one. In the application to a real as- 
trophysical problem this solution, even in the absence of 
rotation, is highly idealized. For example, the presence 
of matter, e.g. in the form of accretion disk, distorts the 
metric. If a static distribution of matter is localized out- 
side the black hole horizon, the spacetime in the vicinity 
of the horizon remains vacuum. We call such a solution 
a distorted black hole. The metric near the horizon of 
a general (not necessary axisymmetric) static distorted 
black hole was studied in 0. 

If the distribution of matter outside a black hole is 
axisymmetric, the metric of the distorted black hole al- 
lows detailed description. The reason is that the vacuum 
metric outside the matter is the Weyl solution. This 
metric contains two functions of two variables. One of 
this function, which has the meaning of the gravita- 
tional potential, obeys the linear Laplace equation in 
a flat 3D space, while the other can be obtained from 
it by a simple integration. Axially symmetric distorted 
black holes were studied in several publications (see e.g. 
0, i, H, i, 0, B 111)- Such axially symmetric distorted 
black holes arise naturally in the models were one of the 
(large) spatial dimensions is compactified. For the gen- 
eral discussion of such solutions in higher dimensions see, 
e.g. [10, 11]. In 4D such caged black hole solution is again 
the Weyl metric. The properties of 4D caged black holes 
were studied in [Hill. 

In the previous studies of distorted black holes the at- 
tention has mainly been focused on the properties of the 
black hole exterior. But any distribution of matter in 
the black hole exterior region distorts the metric not only 
outside the black hole, but also in its interior. The pur- 
pose of this paper is to study this effect. Namely, we 
consider the interior of an axially symmetric distorted 
black hole. In particular, we study the structure of the 
spacetime in the vicinity of the black hole singularity. 
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The paper is organized as follows. In Section 2 we 
collect the equations for the vacuum axisymmetric dis- 
torted black hole in the exterior and interior regions. 
In Section 3 we obtain solution for the interior of dis- 
torted black hole and discuss its properties. An asymp- 
totic form of this solution near the black hole horizon and 
singularity is obtained in Section 4 and Section 5, respec- 
tively. Special examples of exact interior solutions and 
their properties are considered in Section 6. In Section 7 
we consider properties of the interior and singularity of 
4D caged black hole. Section 8 contains summary and 
discussions of the results obtained . Additional technical 
details and calculations used in the main part of the pa- 
per are collected in Appendices. In this paper we use the 
units where G — c = 1, and the sign conventions adopted 
in El. 



II. METRIC OF A DISTORTED BLACK HOLE 

In the absence of distortion a static vacuum black hole 
is described by the Schwarzschild metric 



ds^ = -(1 - 2m/r)dt^ 



dr"^ 



1 — 2m/r 



r^duj^ . 



(1) 



where m is the black hole mass, and dto^ = dO^+siv? 
is the metric on a unit round sphere. In what follows, we 
shall use two other forms of this metric 



ds^ 
ds^ 



cosh tjj ~ 1 
cosh ijj + 1 



dt^ + r7i2(cosh'0 + 1)2 (d^A^ + duj'^ 



^-^2^dt^ + m\l + cosV')'(-d^' + rf^') . (3) 
1 + cos ll) 



The metric ^ is valid outside the Schwarzschild black 
hole horizon, r > 2m, and r = m(cosh'(/; +1), where 
■0 > 0, while the metric ([3]) covers the Schwarzschild 
black hole interior < r < 2m, and r — m(l -I- cos?/'), 
where < 7p < tt. These two metrics are connected by 
the analytical continuation 



tj) — > iip. 



(4) 
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The Carter-Penrose diagram for the interior metric in the 
coordinates {tp, 9) is shown on Figurc[TJ The hnes tp±9 — 
const are null rays propagating within the 2D section 
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{t, (f) = const. It should be stressed that this diagram is 
different from the usual Carter-Penrose diagram for the 
radial sector (i, r) of the Schwarzschild black hole. 




FIG. 1: The Carter-Penrose diagram for the (i/)-0)-sector of 
the Schwarzschild black hole interior. 

Following [3, Q we present the metric for the vacuum 
axisymmetric static distorted black hole in the form 



cosh ij) + \ 
e^^ (de"^ + d^/ 



ds' 



1 — cos V' 



,2U 



1 



COS if: 

2V f ^q2 



dt^ + m^(l + cosV')^e 



-2(7 



sm 



(5) 



(6) 



The metric (O, where U and V are functions of (t/i, 9), is 
valid in the black hole exterior. The metric © , where U 
and V are functions of (V', 8), describes the interior of the 
distorted black hole. The black hole horizon H is defined 
by the equation 

^ = = , < < TT . (7) 

For a regular black hole the functions C/ and V must be 
smooth and finite at the horizon. In particular, these 
functions, which define the geometry of the 2D horizon 
surface, must be continuos across the horizon. The met- 
rics ^ and ^ are regular along the axis of symmetry 
(no conical singularities), provided 



Denote 



D^U = + cot 7p U^^ , 
DeU = U ee + cotetj g. 



(8) 

(9) 
(10) 



Then, the vacuum Einstein equations for the metric ^ 
reduce to the following relations 



(11) 
(12) 

(13) 



Here 

i^W(^,0) 



N 



2 sin ij} cos i> sin Ug 

2 sin -tp cos 9 U^^p — 2 cos -0 sin 6 U^t 



(14) 



N = 



2 sin^ ^ cos 9 UfiU^^ 
— sin ip cos Tp sin 9 {u% 
+ 2 sin ij} cos 9lJ,g — 2 cos ■0 sin 9 U^^ 
sin0(sin^ -0 - sin^ 9)~^ . 



(15) 



(16) 



The functions F'^^ and F'^^^ obey the relation 



F 



{9) 



F 



(17) 



Thus, after solving equation (fTT|) . one can obtain V by 
simple integration 



V{iP,9) 



[F^^)di; + F'^'^U9], (18) 



where the integral is taken along any path connecting 
(Vo,0o) and (V,0). 

The factor N is singular along the lines ijj = 9 and 
ip = TT — 9. Nevertheless, as we shall demonstrate in the 
next Section, the solutions for V for the distorted black 
holes which are regular at the horizon remain smooth and 
regular along these lines. 

Let us focus on equation pT|) . Since D^U — D^^U, 
a solution to this equation can be presented as a sum of 
two solutions, one being odd and the other being even 
function of ip- Because of the presence of factor cot ip in 
D^, this operator is singular at "0 = 0, hence, the regular 
at the horizon solution must be an even function of ip. 
This solution remains real after the Wick's rotation 

Studying the exterior solution 



DgU = 



(19) 



Geroch and Hartle Q demonstrated that if [/ is a regular 
smooth function in any small open neighborhood of H 
(including H itself) which takes the same values, uq, on 
the both ends of the segment H , 



(7(0,0) = C/(0,7r) = Mo, 



(20) 



then the solution is regular at the horizon and describes 
a distorted black hole. 

The surface gravity of the distorted black hole is 



g2uo 

4m 



(21) 



If the distortion source obeys the strong energy condition 
uq has to be non-positive [7]. 
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Let us emphasize that equation ((T9)) for the 'gravita- 
tional potential' U in the exterior region is elliptic, while 
the interior equation (jlip is of the hyperbolic type. This 
is in accordance with the general property of black holes. 
Namely, direction to the singularity in the inner region 
is the direction to the future, and the evolution of the 
metric in this region obeys dynamical equations. After 
solving the equations in the exterior region we obtain 
boundary conditions at the horizon for the inner dynam- 
ical equations. In our particular case the solutions of 
the exterior and interior problems are connected by the 
analytical continuation 

It is convenient to consider dimensionless form of the 
metric dS^, connected to the metric ds^ as follows 



ds' 



Am e 



dS' 



(22) 



Introducing the quantities 

T^Kat, U^U-uo, (23) 
one can write the metric dS^ in the form 

1 -I- cos ^ 4 ' 

X \e^^ (dO^ - dil;^) + sin^ ddcf)^] . (24) 



2^ ^d0^ - dyj^) + sin2 ^-^-^2 



In what follows, we shall study the metric ((24|) and its 
properties. In order to obtain the corresponding charac- 
teristics of the 'physical' solution © it is sufficient to use 
the scaling transformations ((23)) . 



III. INTERIOR SOLUTION 

A. Gravitational potential in the inner region 

Our goal is to study the interior of a distorted black 
hole. To find the metric inside the distorted black hole 
we start with equation (|lip . This equation allows a sep- 
aration of variables U = R{ip)S{9) 



S^ee + cot 0S,e + XS^O, 

R 



(25) 
(26) 



Since the polar points = 0, tt are regular, the functions 
S must be finite at these points. The solutions of this 
eigenvalue problem are 

S' = F„(cos6l), \ = n{n + l), n = 0,l,.... (27) 

Expanding U over the complete set of Legendre polyno- 
mials of the first kind, P„(cos6'), one has 



(28) 



n>0 



Since at the horizon surface tj is finite and regular, one 
must omit infinitely growing at V' = solutions of (j26p . 



Thus, we obtain the following solution for U 

U{iP,e) = ^a„P„(cos^)P„(cos0), (29) 



rs>0 



where a„ are the coefficients called multipole moments. 
For a given value of U on the horizon these coefficients 
are 



a„ = (n + l/2)/ d0 sin 6* C/(0, 6')F„ (cos 61). (30) 
The condition pp]) implies 

a2fc+i = , Q2fc = mq ■ 



(31) 



fc>0 



fc>0 



Since the Legendre polynomials have the symmetry 
property 

P„(-z) = (-l)"P„(z), (32) 
the function U is invariant under the transformation 

/: (V,f?) ^ (7r-V',7r-0), (33) 

that is 

u{-K-ij,Ti-e)^ij{ip,e). (34) 

This relation implies, in particular, that the value of U 
at the singularity ij: = n \s determined by its values on 
the horizon il) — Q 



U{n,TT- 9) = [/(0,( 



(35) 



In other words, there exists an interesting duality be- 
tween the horizon and singularity. It should be empha- 
sized that the functions P^") and F'^'f'\ (see ((HI) - ((T5)) ). 
contain both, symmetric and antisymmetric parts with 
respect to the reflection p3p . This means that the func- 
tion V does not possesses the symmetry (IMl) . Neverthe- 
less, since the function U and boundary conditions ((5|) 
determine V uniquely, the relation (|34p simplifies greatly 
the study of the spacetime structure near the singularity. 
We return to this point in Section V. 



B. Boundary values of U and V 

Let us denote 

u±iO) = ^(±l)"a„P„(cos0) - uo. 

n>0 

It is easy to check that 



u±{0) = u±{-k) = 0. 



(36) 



(37) 
(38) 



4 




71 



FIG. 2: The boundary value of U for the interior of a distorted 
black hole. The dashed line, V = 0, is the horizon, and the 
dotted line, = tt, is the singularity. 



of a free fall of a test particle from the horizon to the 
singularity along the symmetry axis. Consider a particle 
with zero energy and angular momentum which is moving 
along T — const and = 0. The proper time t of its 
free fall from the horizon to the singularity along the 
symmetry axis calculated for the metric (|24|) is given by 



1 



(42) 



The ± signs are for = and 9 = n axes, respectively. To 
illustrate how the distortion of the black hole affects this 
time we consider two simple cases. As the first example, 
we consider the quadrupole distortion when only ao and 
a2 do not vanish. Taking into account that = uq — ai 
(see (|5T|)) one has 



One has (see Figure [2]) 

U{Q,B) = w+(6l), U{7r,e) = u^{9), (39) 
Ui^P^O) - u+ii;), W(^,7r) =w-(V). (40) 

We call these data the boundary value of U. The regu- 
larity of u± follows from the regularity of the distortion 
U on the horizon and the condition ([57)1 . 

The boundary values of V can be obtained from (fT8|l . 
Namely, integrating F^^' along the lines ip = 0, Oq = 0, 
and = n, 9o = one obtains (see Appendix B) 

V{O,0) ^2u+{e), , V{t:,0) ^ ~2u-{e) . (41) 

The boundary values of V are summarized on Figure [31 




FIG. 3: The boundary values of V for the interior of a dis- 
torted black hole. The dashed line, i/; = 0, is the horizon, and 
the dotted line, ip — -k, is the singularity. 



u± = --jCL2 sin V ■ 
The integral ([42[) can be calculated exactly 



(43) 



(44) 



where Iq[z) is the modified Bessel function. A plot of t 
as a function of the quadrupole moment a-i is shown on 
Figure [H (a) . 




FIG. 4: (a): Quadrupole distortion: The proper time, r, 
as a function of the quadrupole moment 02. (b) Octupole 
distortion: The proper time, t+, as a function of the octupole 
moment 03. A plot for r_ can be obtained from this one by 
the reflection 03 —03. 



As a second example we consider the octupole distor 
tion when only a\ and 03 do not vanish. Because of (I3T 
one has a\ = —03, and 



C. Proper time of free fall to the singularity along 
the symmetry axis 

Suppose a function u+(6') for a distorted black hole is 
given. Then the boundary data for U and V described 
above allow one to calculate, for example, a proper time 



u± 



=p— 03 cosi/) sin ip . 



(45) 



To obtain t± we used numerical integration. Plot of 
as a function of the octupole moment 03 is shown on 
Figure [H (b) . The minimal value of r_|_ corresponds to 
a3 ~ —0.9166. The similar plot for r_ can be obtained 
by the reflection aj, — > —03. 
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IV. NEAR HORIZON GEOMETRY 
A. Shape of a distorted horizon 

The form of the horizon surface for the metric p4p is 
determined by the foUowing hnc clement (see also ^]) 

1 — 

This metric is obtained from (|24p as the limit tfj ^ oi 
the metric on 2D section ^ = const of T = const. The 
horizon area is A — At:. 

The Gaussian curvature of the metric du\ is Kj^. — 
R/2, where R is the Ricci scalar curvature. It is given by 
the following expression 

K+ - e-2"+ [1 + (1 - [u'l - 2{u'+ f] - 4zu'_^] , 
= e-2«+ (1 + u+,ge + 3 cot Ou+.g - 2u\ g) , (47) 

where the prime denotes a derivative with respect to z = 
cos 9. 

As the special examples, we consider the quadrupole 
and octupole distortions with functions m+ given by (|43p 
and (I45p respectively. The Gaussian curvature for these 
distortions is 

^ e^^^^i"'" [1 + 802(1 -5 cos^e) 

- ISa^ cos^ 6* sin^ 6*] , (48) 

= igS^^cosesin^e lOaa cos 6'(9cos2 61 -5) 

- 25a^sin2 6i(l -Scos^e)^] . (49) 

For the quadrupole distortion the Gaussian curvature be- 
comes negative at both of the poles, 9 = and 9 = tt, 
for 02 > 1/12. Similarly, for the octupole distortion, the 
Gaussian curvature becomes negative at one of the poles 
for joal > 1/20. It means that for these values of the pa- 
rameters the horizon surface of the distorted black hole 
cannot be isometrically embedded in a flat 3D space (see 
e.g. [1^ and references therein). For 02 < 1/12 (in the 
quadrupole case) and jaal < 1/20 (in the octupole case) 
isometric embeddings are possible. 

To construct the embedding we consider a surface 

p = p{e) , z = z{e) (50) 

in 3D Euclidean space with the metric 

df = dz^ + dp^ + p^dcf)"^. (51) 
The geometry induced on this surface, 

dP = {z% + p^g)d9^ + p'df , (52) 
coincides with the horizon surface geometry (j46p if 

p = e-"+sin6l, z= d9Q, (53) 
Je 

q2 ^ g2«+ _ g-2«+ (^Qg 9-u+^g sin 9f. (54) 




FIG. 5: The shape of the horizon surface of the distorted black 
hole. The embedding diagram for the horizon surface can be 
obtained by rotation of the curves on the plots around the ver- 
tical axis. The left plot (a) shows the rotation curves for the 
quadrupole distortion of 02 = 1/12 (line 1), and a2 = —1/12 
(line 2). The right plot (b) shows the rotation curves for 
the octupole distortion of a\ = —as = 1/20 (line 1), and 
ai = —az = 1/6 (line 2). The region of negative Gaussian 
curvature is schematically illustrated on plot b by the dashed 
horizontal line. The rotation curves for positive octupole mo- 
ments 03 can be obtained by reflection of the lines on plot b 
with respect to the horizontal axis. Dashed lines on both the 
plots are round circles of radius 1. 



Figure [5] shows the embedding diagrams of the distorted 
event horizon surface for the quadrupole and octupole 
distortions. 



B. Kretchmann invariant 

The bmction u+(0), which specifies the geometry of 
the horizon surface, uniquely determines the geometry 
of the black hole interior. In particular, one can obtain 
expansion of U and V at the vicinity of the horizon (see 
Appendix B). The first two terms of this expansion in 
the powers of "0 are 

U = + ■ • ■ , (55) 

V = 2u+- i('u^^' - u\ g + 2cot6lw+,e)V'^ + . . . . 

(56) 

Here and later we use the dots . . . for the omitted terms 
of higher order in -0. We also defined 

u^l\9) = ^(±l)"a„n(n+ l)P„(cos6'). (57) 

n>0 

In this approximation the metric near the black hole hori- 
zon reads 

dSl = A+dT^ + B+{d9^ - di}^) + C+d(f>^ , (58) 
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2u\ 
1 



1)^ 



(2) 



4 cot ^u+^i 



(59) 



C+ = ^e"2"+sin2 6i[2 + (m 



(2) 



This expansion, for example, can be used to determine 
the value of the Kretchmann scalar K, = Rap-ysR"^^^ at 
the horizon surface. Using ([55)1 calculations give 



-4n 4 



'l + u+,ei0 + 3 cot 6u+.e — 2u\ 



(60) 



For undistorted Schwarzschild black hole /Csd, = 12. 
Figure [6] illustrates the ratio, k — /C+//Csch_+, of the 
Kretchmann scalars of the distorted and Schwarzschild 
black holes. 




FIG. 6: The ratio k = K.+ /K,scii j_, of the Kretchmann scalar 
>C+ on the horizon of the distorted black hole to its undis- 
torted value k = A^sch + . Curves on plot (a) show k for the 
quadrupole distortion of a2 = 1/12 (line 1), and a2 = —1/12 
(line 2). Similar curves on plot (b) show k for the octupole 
distortion of ai = —03 = 1/20 (line 1), and ai = —03 = 1/6 
(line 2). The dashed horizontal lines at fc = 1 correspond to 
the Schwarzschild black hole. 



obtained from its asymptotic expansion near the horizon 
-0 = 0. After this, the relation (|18p allows one to find the 
expansion of V . The expansions for U and V near the 
singularity are given in Appendix B. Using these expan- 
sions one can obtain the asymptotic form of the metric 
p4)) at V- = 71- - V' ^ 



(62) 



B_ = 



96 



-[6+ (3u_,ee + 3cot6'u_,e - 1) V'? 



-V'l[6- (9u_,fl 



(63) 



3 cot eu^^g - 6(u_,e)2 + l)tp'i + 



-2u- 



■ sin^ dil;'^ 



96 



Expressions ([63|l are sufficient to calculate the Kretch- 
mann scalar near the singularity up to the second order 
in ip^ corrections 



r(2) 



49152 ei2"- / 1 „ ,^f2) 



(64) 



1 + 3u^^gg — 6{u^,g)'^ — 3 COt 9u^,g . (65) 



Higher order terms can be obtained by using the relations 
given in Appendix B. In the absence of distortion, when 
u_ = 0, the Kretchmann scalar does not depend on 9 



/Csoh 



49152 



(66) 



This is the value of /C_ for the Schwarzschild geometry. 
It can be shown that the metric of a distorted black hole 
is of type D near the singularity (l^ • 



Comparing this result with (|47p one can arrive to the 
following relation valid on the horizon surface of a dis- 
torted black hole 

1C+ = 12Kl . (61) 

This relation is valid not only for axisymmetric case but 
also for an arbitrary static distorted black hole (see Ap- 
pendix It is possible to show that the metric of the 
distorted spherical black holes is of type D both on the 
horizon and on the axis of azimuthal symmetry, . 

V. SPACETIME NEAR THE SINGULARITY OF 
A DISTORTED BLACK HOLE 

A. Asymptotic form of the metric 

Because of the symmetry property the asymp- 

totic form of U near the singularity ijj = n can be easily 



B. Stretched singularity 

For the Schwarzschild geometry the metric near the 
singularity 

dSl^-^i;td^^_ + ^dT' + ^du;\ (67) 

can be written in the form 

dSl « -dr^ + (^dT^ + ^^^d.^ • (68) 

Here r = —ip^/12 is the proper time of a free fall to 
the singularity along the geodesic T,9,(f> = const. The 
quantity r is negative, and it reaches at the singular- 
ity. The metric ([68)) has the Kasner-like behavior with 
indices (—1/3,2/3,2/3). It describes a metric of collaps- 
ing anisotropic universe, that shrinks in 9-(f> directions, 
and expands in T— direction. 
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The Kretchmann invariant as a function of the proper 
time has the foUowing asymptotic form 



/Co 



64 
27t4 



(69) 



This relation shows that the surface of constant A^sch _ is 
at the same time a surface of constant r. 

Spacetime in the region where the curvature is of or- 
der of the Planckian curvature requires quantum gravity 
for its description. For the Schwarzschild geometry at 
the surface where ICs^h _ ^ Ipf the proper time r is of 
order of the Planckian time Tpi. Since one cannot rely 
on the classical description in this domain, it is natu- 
ral to cut the region where the curvature is higher than 
the Planckian one and to consider its boundary as the 
stretched or 'physical singularity '. For the Schwarzschild 
metric the stretched singularity surface has the topology 
X S^. Its metric is a direct sum of the metric of a 
round two-sphere and a line. 

What happens to the stretched singularity when the 
metric of the black hole is distorted? To answer this 
question we use the asymptotic form of the metric near 
the singularity, ([5^ . Consider a timelike geodesic lying 
on the 'plane' T = const, (j) = const. We call such a 
geodesic 'radial'. It can be shown (see Appendix C) that 
a 'radial' geodesic is uniquely determined by the limiting 
value d of its angular parameter 9 at which it crosses 
the singularity. Denote by t the proper time along the 
'radial' geodesic to its end point at the singularity. In 
coordinates (r, d) the metric dS^ is given by ([55]) where 
dcj^ is replaced by 



dO^ 



(70) 



We can use (r, d) as new coordinates in the vicinity of the 
singularity. Relations (|C14p . (|C15|) connect these 'new' 
coordinates with the 'old' ones {ip-,9). The Kretchmann 
scalar (|64|. ((65l) in 'new' coordinates reads 



64 
27^4 



[l + ^L'V2/3 + 0(r4/3)], (71) 



^L'^ = i(12)2/3e2"-W[l + 3?/_,^^ 
— 6(u_^^)^ — 3cot0w-,ij]. 



(72) 



The expansion (|7ip coincides in the leading order with 
(|69)) . Hence, in the presence of distortion surfaces of 
equal )C are again (in the leading order) surfaces of con- 
stant T. 



C. Shape of equi-curvature surfaces 

A surface S- where the Kretchmann scalar has con- 
stant value /C_ = JCq- In the vicinity of the singularity 
(in the leading order in ?/;_) tp^ and on I]_ arc related 
as follows 



K- = (49152//Co) 



1/12 



(73) 



and one has the relation tp-^e = ''P-U-fi. Consider the 
induced geometry on Using this relation one can 

conclude that dip'^ term in (|62p gives quadratic in ■;/)_ 
corrections only. Neglecting all such terms in ([5^ we ob- 
tain the following expression for the leading asymptotic 
for the induced metric dP on S_ 



dll « ^dT'^ + —da'i 
nt 16 



(74) 



where dcr^ is given by ([70]) . The surface S- has the same 
topology X S'^ as in the absence of distortion, but its 
geometry is different. This difference manifests itself in 
the shape of T = const 2D surfaces. The information 
about the shape is encoded in the 2D metric rfcr^. The 
total area of the surface is 47r. The metric ([70)1 can be 
obtained from the horizon metric da\, (j46[) . by a simple 
change Uj^ — > — u_. Under this transformation the even 
multipole coefficients are invariant, while the odd coef- 
ficients change their sign (see (IM)) '). In particular, the 
embedding diagrams for the metric (j70p are the same as 
shown on Figure [5] (a) , with the same value of 02 , and 
with the opposite value of the octupole moment 03, as 
shown on Figure [S] (b). 



VI. EXACT SOLUTIONS 

For any given set of multipole coefficients that de- 
termine J7, the corresponding function V can be found 
explicitly in terms of elementary functions. Since the 
general expression for V is rather cumbersome we do not 
present it here. Instead, we consider special case of the 
quadrupole and octupole distortions when the gravita- 
tional potential U is of the form 

3 5 
U = --02^3/2 - -03 cos V cos 61^5/2, (75) 



where 



Vq — siv? ip + sin^ 6 — q sin^ ip sin^ I 



(76) 



Substituting this solution into equations (|T4l) . (1151) . ^^d 
integrating according to equation (jlSp we derive 



V — sin^ 6'[— 3a2 cos ■(/; -I- -03 cos 0(1 — 3 cos^ V')] 
^ ,2 /)r„2-i 



+ - sin^ V sin^ f [02^22 + 20203 V23 + V33] , (77) 



where 



V22 = 9[l-rg/s], 

V23 = 15cosV'cos6i[l - 3/2P3/2] , 

V33 = ^[4-127'39/24 - 187'2' + 27Pi\/6] 



(78) 



Exterior metric for a black hole distorted by a quadrupole 
field was derived in ^171] . 
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Using GRtensorll package we calculated the Kretch- 
mann scalar K. for this distortion. We made the calcula- 
tions for both the inner and external regions. It is easy to 
check that IC^s in these regions are related by the analyt- 
ical continuation ([4]). Figures [7] and [8] show the contour 
lines of K, for the quadrupole and octupole distortions re- 
spectively. In order to plot both the exterior and interior 
regions simultaneously we introduce new coordinate 

cos tp , for £ (0, tt) ; 

cosh tp , for ip (0, oo) . 

The sector y G (—1,1), z £ [—1,1], where z — cosd, 
covers the inner region, and the sector y G (l,oo), z G 
[—1,1] covers the exterior of the black hole. 



-0.5- 
0.0- 
y 0.5- 
1.0- 
1.5- 




0.0- 
y 0.5 ; 
1.0^ 
1.5- 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 

Z Z 

a b 

FIG. 7: The contour lines of K, for the quadrupole distortion 
of 02 = 1/12 (a), and 02 = —1/12 (b). The horizontal line 
y = 1 represents the event horizon. 




FIG. 8: The contour lines of IC for the octupole distortion 
of ai = -as = 1/20 (a), and ai = -as = 1/6 (b). The 
horizontal line y = 1 represents the event horizon. 



VII. INTERIOR OF A CAGED BLACK HOLE 

In this Section we apply the obtained results to a spe- 
cial case of the so called 'caged' black hole. Such a black 



hole is a solution of the vacuum Einstein equations for 
a spacetime where one (or more) spatial dimensions are 
compactified. Caged black holes in 4D spacetime were 
discussed in fTs*] (see also [l^, As the result of the 

compactification the event horizon of the black hole is 
distorted. The metric is axisymmetric and is a special 
case of the Weyl solution. The value of U at the horizon 
found in [l^ [ equation (62)] gives 



t/(0, 9) = ^ ln(47r) + - In 



/ 



^J. + z 



f 



(79) 

where jz] < ^ and < ^ < tt. In this case z — /icos6', 
/i is a dimensionless parameter equal to the ratio of the 
black hole mass, m, to the radius of the compactification, 
L, = m/L, and 



f{x) = —X sinx r^(a;/7r). 



(80) 



The function f{x) has the following properties: 
/(0) = 1, /(vr/2) = l, /(^)=0, 



n) 



and at the interval < a; < tt it can be approximated by 
a linear function 



/(^) 



1 



X 
TT 



(82) 



with accuracy of 1% . The coefficients a„ for this solu- 
tion U can be obtained from ([50]) . Let us emphasize that 
since the function U{0, 9) is invariant under the transfor- 
mation 9 ^ TT — 9, the coefficients a„ for odd n vanish 
(see equations ([BG)) . ([BT)) ). This implies that u_ = u+ 
and the boundary value of U at ip ^ n coincides with the 
boundary value of this function at the horizon, -0 = 
(see equation ([39])), and we have (see equation (pO|) ) 



Uo 



^ ln(47r) 

TT 



I ln[/(A^)]. 



From equations (|79|) and (|83p we derive 

1 / 47r(7r - /i) + ^2 sin^ 1 
u{9) = u±{9) « - In f 



2 V 47r(7r - fi) 
The metric P5|) on the surface of the horizon is 



(83) 



(84) 



dai ^ 1 



1 



^2 sin^ 9 
47r(7r — fi) 

^2 sin^ 9 
47r(7r — fi) 



d9^ 



sm^9d(l)^. (85) 



Using equations ((60|) . ((6T|) . and ((7T|) . (jTS]) we can calcu- 
late the Kretchmann scalar at the horizon surface of a 
caged black hole 



/C+ w 12 (47r(7r - ^i)) 



4 (47r(7r - ^) + ^2 _^ 3/i2 cos^ Of 



(47r(7r - fi'^ sin^ 9) 



(86) 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 



FIG. 9: (a) The shape of the distorted event horizon sur- 
face for the caged black hole. (1): fi = 2/37r, (2): jj, — G/Yn. 
The dashed circle corresponds to the Schwarzschild black hole. 
The embedding surface is obtained by rotation of these curves 
around the vertical axis, (b) The Kretchmann scalars ra- 
tio k = /C+//Csch + on the horizon for the same values of fi 
(lines 1 and 2) . The dashed horizontal line corresponds to the 
Schwarzschild black hole. 




FIG. 10; The shape of the distorted 'physical singularity' for 
the caged black hole of fi — 2/37r (line 1), and fi = 6/7n (line 
2). The region of negative Gaussian curvature is schematically 
illustrated by the dashed horizontal lines. The dashed circle 
corresponds to the Schwarzschild black hole. The embedding 
surface is obtained by rotation of these curves around the 
vertical axis. 



and in the vicinity of its 'physical singularity ' 
64 



.[1 + ^(2)^2/3^0(^4/3)] 



(2) 



27t4' 



(87) 



-I- 3Ai'*sin^6i(13sin2 6l - 16)] 

X [I67r(7r - fi) (47r(7r - fi) + fi^ sin^ 9)] .(88) 



Applying equations ([53|) . ([54|) . and the transforma- 
tion u+ — > — u_ respectively we can construct embed- 
ding of these surfaces. Figure [9] illustrates the shapes of 
the distorted event horizon surface and the ratio of the 
Kretchmann scalars , k = 1C+ /A^sch + > of the caged and 
Schwarzschild black holes. The metric ([70|l on the surface 



of the stretched singularity is 

/i2 sin^ e 



dot 



1 



1 



47r(7r — /i) 
sin^ 6 
47r(7r — /i) 



dO' 



sm 



(89) 



The shape of distorted 'physical singularity ' is illustrated 
on Figure [TOl 



VIII. DISCUSSIONS 

Let us summarize the obtained results. We considered 
the geometry of static vacuum axisymmetric distorted 
black holes. We focused mainly on the properties of the 
horizon and interior of such black holes. The geometry 
of a distorted black hole is uniquely determined by the 
'gravitational potential' U which is a solution of the 3D 
flat Laplace (in the exterior) or d'Alembert (in the inte- 
rior) equation. After solving this 'master' equation, the 
second function V , which enters the metric, can be ob- 
tained by a simple integration. 

The 'gravitational potential' U can be written as a su- 
perposition of the Schwarzschild potential J/sch and the 
distortion potential U . The distortion U is determined 
by the values of the multipole coefficients obeying the 
constraints ([?T|) . The distortion potential in the black 
hole interior possesses a remarkable discrete symmetry 
([M)) which relates the value of U in the vicinity of the 
singularity to its value at the horizon. Thus, the func- 
tions u±{9), (see dH])), determine both, the shape of the 
horizon and the leading asymptotics of the metric and 
curvature invariants near the singularity. 

Qualitatively, the shape of the event horizon surface of 
distorted black hole is similar to the shape of equipo- 
tential surfaces in the linearized (Newtonian) gravity. 
Namely, consider a point-like mass M . In the presence 
of a quadrupole distortion its Newtonian gravitational 
potential reads (up to constant oq) 



$ = h A$ . 

R 



2 ^ ' 



(90) 



where R = \/ + y'^ + is the radial distance from the 
mass M , and a2 is the value of the quadrupole moment. 
For positive a2 such a distortion is generated, for exam- 
ple, by a ring of mass m and radius d ^ R located in the 
equatorial plane. For such a ring a2 = m/(2d^). Sim- 
ilarly, a negative 02 is generated, for example, by two 
point masses m located on the axis of symmetry on the 
opposite sides of the mass M at the distance d ^ R. In 
this case a2 — —2m/d^. We consider R ^ M and as- 
sume that the distortion A$ is small. Then the change 
5R in the position of the equipotential surface for (^0)) 
with respect to the position of the unperturbed surface 
of i?o = const is 



5R 



0.2 Rq 

' 2M 



(3 008^6'- 1). 



(91) 
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Thus the quadrupole distortion deforms the equipotential 
surfaces and makes them either oblate (for 02 > 0), or 
prolate (for 02 < 0). This property is similar to the 
property of the horizon surface for the distorted black 
hole (see Figured (a)). 

It should be emphasized that the linear approximation 
is not sufficient for the 'explanation' of the Kretchmann 
invariant properties. Really, in the linear approximation 



-(1 + 2$)dt2 + (1 - 2^){dx^ + dy^ + dz^) , (92) 



the Kretchmann scalar is 



M{M + 2m) 

fl6 



Its variation under the small distortion A$ is 

2 Af ^ 
5fC = -^5R. 



(93) 



(94) 







Hence in the weak field approximation, /C is larger at the 
points where SR < 0, such as at the pole of the oblate 
equipotential surface (for 02 > 0), and in the equatorial 
points of the prolate surface (for 02 < 0). This behavior 
of JC in the weak field limit is opposite to the behavior 
of AC on the horizon of the distorted black hole (see e.g. 
(ICTj) ). This difference demonstrate that non-linear effects 
and the spatial curvature are important near the horizon. 

The property (j6ip has an important consequence for 
caged black holes discussed in Section IVIII For /i close 
to TT, when the 'north' and 'south' poles of the caged 
black hole are close to one another, the Gaussian cur- 
vature (and hence the Kretchmann invariant) becomes 
large at the poles. In other words, in the infinitely slow 
merger transition, the region of a very high curvature 
'leaks' through the horizon in the vicinity of the black 
hole poles. When this curvature reaches the Planckian 
value, one can say that the 'physical singularity' (as de- 
fined in subsection B of Section V) becomes naked. This 
may indicate that during the phase transition between 
black-hole and black-string phases one can expect a for- 
mation of a naked 'physical singularity'. Whether this 
conclusion remains valid for higher dimensional caged 
black holes and beyond the adiabatic approximation is 
an interesting open question. 



Acknowledgments 

This research was supported by the Natural Sciences 
and Engineering Research Council of Canada and by the 
Killam Trust. 



static distorted black hole is related to the Gaussian cur- 
vature of the 2D horizon surface Kq as follows 



(Al) 



Let us emphasize that this relation is valid for an arbi- 
trary (not necessarily axisymmetric) distorted black hole. 
To establish this relation we use the results of paper [2]. 

The metric near the horizon of a static vacuum dis- 
torted black hole can be written in the Israel coordinates 
ni as follows 



ds^ = -Xdt^ 



dX^ 



habde^de^ 



(A2) 



where —X — is the square of the timclike Killing vec- 
tor it = dt, 

K = ipi/z ^ -2^<j|f,^"l'' = = X-\VXf . (A3) 

Here (. . .)a\b means the covariant derivative with respect 
to 2D metric hab, V and □ are the operators associated 
with this metric. Using the properties of the Killing vec- 
tor one can show that in the Ricci flat spacetime the 
following relations are valid 



dxhab = K ^kab, dxK = -l-k 



'^a b ' 



k — h kab : 



XdxK. = -k 



1 



1 



(A4) 
(A5) 
(A6) 



tkS^^ - -(k-^). 1° 
2 " 2^ ^1° 

{AKr^Xiik.dk^" - e)6^ + 2fcC] ,(A7) 



K = Kk- ]^X{kabk'''' 



(A8) 



Here K is the Gaussian curvature of 2D surface X = 
const. The Kretchmann invariant for the metric (|A2p 
can be written as follows (see equation (4.18) of 0) 



IC = SK^ikabk'"'' + k^ + 2(Xk2)-1k|^kI"] 



(A9) 



One can show that quantities hab, kab and k are finite at 
the horizon and its vicinity they allow a representation 
in the form of the Taylor series 

K = Kq + KiX + 0{X'^) , Ko|a=0, (AlO) 



C = lko6a' + 0{X), 
K = Ko + 0{X) , Ko = Kofco . 



(All) 
(A12) 



Equations (jAlOp imply that the third term in the squared 
brackets of (jA9|) vanish at the horizon X = 0. Simplify- 
ing the other two terms by using (jAlip and (jA12p one 
obtains (Ml. 



APPENDIX A: KRETCHMANN INVARIANT 
AND GAUSSIAN CURVATURE ON THE 
HORIZON OF A STATIC BLACK HOLE 



APPENDIX B: ASYMPTOTIC EXPANSIONS 
NEAR THE HORIZON AND SINGULARITY 



In this appendix we show that the Kretchmann invari- 
ant /C — R^uXpR^'^^'^ calculated at the horizon of a 4D 



The solution (|29|) can be used to find the asymptotic 
behavior of U near the horizon, -0 = 0, and singularity, 
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Tp = IT. To deal with both cases simultaneously, we de- 
note V'+ = "0 s-iid "0- = — ip- The function U is an 
even function of ipa- (c = i), and it has the following 
expansion 



(Bl) 



Here J/i-^"'' are functions of angle 6. The operator in 
([5]) has the same form Dg. for both the variables ip„ 

D, = dl^+coii^,d^^. (B2) 

Using the series expansion for cot ip^ 



cot V'.T = 0g Ml - XI '=2™V'?"] , 
m— 1 

(_l)m-l22m52^ 



C2n 



(2m)! 

where B2m are the Bernoulli numbers 



(B3) 
(B4) 

(B5) 



the relation 

oo 

D^^l- = 4nV^^"-') - 2n X c,^i^T+"^-'\ (B6) 

m— 1 

and equation one obtains 









(B7) 




1 , 


cot 6'uCT^e) , 


(B8) 






, + cot0;7(') + ^[/(2)j , 


(B9) 


rr(2n+2) 
'-^ cr 


1 


-m + l)c2„C/(2("~™+i))] 




4(n+ 1)2 

n 

711—1 





(BIO) 

Similarly, the asymptotic expression for V near the 
horizon and singularity can be written in the form 



(Bll) 



ri=0 



where Va'^^'^ are functions of the angular variable 0. Sub- 
stituting expansion (|Bip into equation (|18p one can de- 
termine the functions vi^""* . The first three of these func- 



tions are 



(B12) 



yj^) = ^{a[24t/^4)+6(l + 2cot2 0)[/^2) 

- 6 cot eufl - (5 + 6 cot^ e) cot e u„^g\ 

- 12cot6»[/PMg,e + 12(C/P)2}. (B14) 



APPENDIX C: GEODESIC MOTION NEAR THE 
SINGULARITY 

For a free particle moving in the black hole interior 
there exist two integrals of motion connected with the 
spacetime symmetry 



(CI) 



The first has the meaning of the conserved momentum 
along T axis, and the second one is the angular momen- 
tum. 

Consider a point (■0° , 0°, T", 0°) near the singularity 
of a distorted black hole. What is the proper time r° 
required to fall from this point to the singularity? This 
time depends on the value of E and L. We consider 
the proper time for the special value of these parameters 
E = L — 0. In this case, for a moving particle T — const 
and L — const. For the Schwarzschild geometry this is 
the radial motion. We also call the motion in the interior 
of a disported black hole 'radial' when E — L — 0. This 
type of motion is a geodesic in 2D metric 



obtained by dimensional reduction from 
Let us denote 

a=i(lni?_),^_, /3=i(lnB_) 
Then the Christoffel symbols for the metric dj'^ are 



-^8 _ -r^- 



The geodesic equation 

d^x'' „ dx" dx^ 
dT^ dr dr 

in the metric (IC2I) takes the form 



V'- + a{ip'i + e^) + 2/3V^_6l = , 



L -1-/3(02 +^2) + 2a0_( 



0. 



(C2) 

(C3) 

(C4) 
(C5) 

(C6) 



(C7) 
(C8) 
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Here the overdot denotes derivative with respect to the 
proper time r. These equations obey the constraint 

B^iijl-e"^) = 1, (C9) 

that is the normahzation condition, Ui^^u'^ = —1, for 4- 
velocity. 

Using (|63p for the metric near the singularity we have 

hiS_ =4hi!A_ -6u_ -\nl6 + 0{i;l). (CIO) 

Thus, in the leading order a « 2/?/'_, and (3 ~ —Su^^g. 

In the leading order the geodesic equations (|C7P - (|C8p 
and the constraint (jC9|l take the form 

^-■ip- + 2{^'i + _ Qu^^g^j^^^e « 0, (Cll) 

il^Je -Zu-^gij-{i?_ + ^2) + 4,/,_^~o, (C12) 

e-s^-V'i (^A^ - ~ 16. (C13) 
According to (jClOp . the order of approximation in the 



geodesic equations corresponds to the order of approxi- 
mation of the metric ((62|) . 

We use the ambiguity in the choice of r to put t = at 
the singularity for each of the 'radial' geodesies approach- 
ing the singularity. Since r grows along geodesies di- 
rected to the singularity, it is negative before the geodesic 
approaches the singularity. The point r = is a singu- 
lar point of the equations (|Clip - (jC13|) . To find approx- 
imate solution to the geodesic equations one can apply 
the method of asymptotic splittings described in [1^. A 
'radial' geodesic approaching the singularity is uniquely 
determined by the limiting value ■& — 8{t = 0). The 
asymptotic expansion of tp^ and 6 near r = is of the 
form 

^- = + + 0(f^/3), (C14) 

= i9+iu_,^(i?)f2/3 + 0(f4/3)^ (C15) 

where f = -12e3"-('')T. 
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